Triaxial Analytical Potential-Density Pairs for 

Galaxies 

Daniel Vogt*and Patricio S. Letelier^ 
Departamento de Matematica Aplicada-IMECC, Universidade Estadual 
de Campinas, 13083-970 Campinas, S. P., Brazil 

February 1, 2008 



Abstract 

We present two triaxial analytical potential-density pairs that can 
be viewed as generalized versions of the axisymmetric Miyamoto and 
Nagai and Satoh galactic models. These potential-density pairs may be 
useful models for galaxies with box-shaped bulges. The resulting mass 
density distributions are everywhere non-negative and free from sin- 
gularities. Also, a few numerically calculated orbits for the Miyamoto 
and Nagai-like triaxial potential are presented. 

1 Introduction 

There are several three-dimensional analytical models in the literature for 
the gravitational field of different types of galaxies and galactic components. 
Jaffe PQ and Hernquist [2] discuss models for spherical galaxies and bulges. 
Three-dimensional models for flat galaxies were obtained by Miyamoto and 
Nagai [3] and Satoh [4j; de Zeeuw and Pfenniger [5] presented a sequence 
of infinite triaxial potential-density pairs relevant for galaxies with massive 
haloes of different shapes. Long and Murali [6] derived simple potential- 
density pairs for a prolate and a triaxial bar by softening a thin needle with 
a spherical potential and a Miyamoto and Nagai potential, respectively. See 
[7] for a discussion on other galactic models. There also exist several general 
relativistic models of disks, e. g., [8]-|15j. A general relativistic version of 
the Miyamoto and Nagai models was studied by Vogt and Letelier |16j . 
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Although axial symmetry is a common approximation to many disk 
galactic models, recent statistics of bulges of disk galaxies (SO-Sd) have 
revealed that nearly half of them are box- and peanut- shaped [T7]. About 
4% of all galaxies with box- and peanut- shaped bulges have so called "Thick 
Boxy Bulges": they are box shaped and are large with respect to the di- 
ameters of their galaxies [18J. In this work we consider two simple triax- 
ial potential-density pairs whose mass density distributions are box-shaped 
along one axis. The first is constructed by applying on each Cartesian coor- 
dinate a transformation similar to that used by Miyamoto and Nagai; this 
is done in subsection 12.11 The second pair is a triaxial generalization of one 
studied by Satoh, and is presented in subsection !2.21 In section [3] some orbits 
for the Miyamoto and Nagai-like triaxial potential are exhibited. Finally, 
the results are discussed in section [H 



2 Triaxial Models for Galaxies 

In the following potential-density pairs, the mass-density distribution is ob- 
tained directly from Poisson equation 

' (*,*» + *, W + *,zz) , (1) 



r A-kG 

where $(x,y,z) is the gravitational potential. 

2.1 Triaxial Miyamoto and Nagai-like Model 1 

We start with the gravitational monopole potential in Cartesian coordinates, 

Gm 



yj x 2 + y 2 + z 2 



and apply the transformations x — > ai + \J x 2 + b 2 , y — > a 2 + \Jy 2 + b\ and 
z — > 03 + y 1 z 2 + 63, where the at, hi are non-negative constants. Using 
equation ([1]), we obtain the following density distribution 

! {^V[(ai + 2 («l + 30 



-01(02 + t]) 2 + oi(l + x) 2 ] + [(«2 + r/) 2 (a 2 + 3r?) + a 2 (ai + 2 

-a 2 (l + x) 2 ] + Hfv 3 [(1 + X) 2 (l + 3x) + (01 + 2 + («2 + r?) 2 ] } , (3) 
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Figure 1: Constant density curves of equation ([3]) on the planes (a) z = 0, 
(b) x = and (c) y = with parameters a\ = 1.0, b\ = 1.0, <i2 = 0.5, 
62 = 1-0 and 63 = 0.5. The contour levels in (b) and (c) are the same as 
shown in (a). 



where the variables and parameters are rescaled in terms of 03: x = x/a^, 

V = y/az, z = z/a 3 , a,i = cii/a 3 , hi = bi/a 3 , p = a\pjm, and £ = J x 2 + b\, 

V = \jy 2 + H and x = \j z 2 + 63. 

The density distribution equation ([3|) is always non-negative and free 
from singularities. For a\ = ai = the potential-density pair is axisymmet- 
ric with respect to the z axis, and in particular if b\ = 62 = we recover the 
Miyamoto & Nagai model 1. Figures QJa)-(c) show some isodensity curves 
of equation ([3|) on the planes (a) z = 0, (b) x = and (c) y = with 
parameters a\ = 1.0, b\ = 1.0, 02 = 0.5, 62 = 1.0 and 63 = 0.5. From a 
top view, the matter distribution is box-like shaped, whereas from a lateral 
view matter seems more flattened in a disk- like manner. Plotting some iso- 
density curves with other values of the parameters allows us conclude that 
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larger values of a\ and a 2 lead to larger deviations from axisymmetry, and 
the degree of flatness with respect to the plane z = depends on 63 as in 
the original Miyamoto and Nagai models. 

2.2 A Triaxial Satoh-like Model 

Satoh [3] derived a family of three-dimensional axisymmetric mass distri- 
butions by flattening the higher order Plummer models of order n. When 
n — > 00 the potential takes a rather simple form 

Gm 



1/2' 



(4) 



x 2 + y 2 + z 2 + a(a + 2Vz 2 + b 2 ) 
We propose the following triaxial generalization of the above potential: 



<£> = —Gm 



+ 01(01 + 2Jx 2 + b\) +y 2 + a 2 (a 2 + 2Jy 2 + b 2 2 ) + 



+a 3 {a 3 + 2J z 2 + b 2 3 ) 



-1/2 



(5) 



The corresponding mass density distribution follows from equation ([T]) 

1 



P 



^ev 3 X 3 [{x 2 + 01(01 + 2^) + y 2 + a 2 (a 2 + 2rj) + z 2 + 1 + 2 X f 2 
{aib\rfx Z [x 2 + y 2 + z 2 + (01 + 2e)(ai + 3^ + a 2 (a 2 + 2j?) + 1 + 2 X ] 
+a 2 62? 3 X 3 [x 2 + y 2 + z 2 + (o 2 + 2r/)(a 2 + 3r?) + 01 (01 + 2^) + 1 + 2 X ] 
+b 2 CW [x 2 + f + z 2 + {l + 2 X )(1 + 3 X ) + 01(01 + 20 + a 2 (a 2 + 2r/)] } , 

(6) 

where the variables and parameters were rescaled as in subsection 12.11 The 
density distribution equation is also non-negative and free from singu- 
larities. For ai = a 2 = we recover the original Satoh model. In figures 
E^a)-(c) we display some isodensity curves of equation on the planes (a) 
z = 0, (b) x = and (c) y = with parameters a\ = 1.0, b\ = 1.0, a 2 = 0.5, 
b 2 = 1.0 and 63 = 0.5. 



3 Orbits in triaxial potential-density pairs 

In this section we briefly discuss some orbits calculated numerically for the 
triaxial Miyamoto and Nagai-like potential-density pair presented in subsec- 
tion 12.11 We first consider orbits on the z = plane. Since the potential 
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Figure 2: Constant density curves of equation © on the planes (a) z = 0, 
(b) x = and (c) y = with parameters ai = 1.0, b\ = 1.0, a2 = 0.5, 
62 = 1-0 and 63 = 0.5. The contour levels in (b) and (c) are the same as 
shown in (a). 
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Figure 3: Radial coordinate R = R/a^ as function of polar coordinate (p for 
motion with triaxial potential (solid curves) and with axisymmetric potential 
(dashed curves). Parameters: a\ = 1.0, b\ = 1.0, 0,2 = 0.5, 62 = 1.0 and 
63 = 0.5. Initial conditions: R = 3.0, R = 0, p = and p was calculated for 
an energy E = —0.15 in (a) and E = —0.10 in (b). 



is a function of the polar angle, there is no conservation of the z compo- 
nent of the angular momentum. In figures El^a)-(b) we compare the orbits 
for the triaxial potential (solid curves) and for the axisymmetric potential 
(dashed curves). The parameters used were the same as in figure [Hand for 
the axisymmetric case a\ = 0,2 = 0. In the equations of motion time was 
rescaled as t = (Gmj 'a\) l ^ 2 t and the conserved mechanical energy per unit 
mass E was rescaled as E = (Gm/as)^ 1 E. Initial conditions in cylindrical 
coordinates were set as R = 3.0, R = 0, (p = and ip was determined for a 
total energy E = —0.15 in figure EUa) and E = —0.10 in figure E^b). In the 
axisymmetric case the orbits shown are precessing ellipses. When axisym- 
metry is destroyed, so is the regular oscillation of the radial coordinate. At 
a higher energy, both orbits are qualitatively almost identical. 

Figures [Ha)-(d) show an example of three-dimensional orbit for (a) ax- 
isymmetric potential and (b) triaxial potential. In figures IHc)-(d) the radial 
coordinate R and the coordinate z, respectively, are plotted as function of 
time t for the triaxial case (solid curves) and axisymmetric case (dashed 
curves). Here the parameters used were the same as in the previous exam- 
ple; the initial conditions were R = 3.0, R = 0, ip = and the imposed 
initial constrait z = 0.5R<p was calculated for an energy of E = —0.15. As 
before, motion with triaxial potential is more irregular compared with the 
axisymmetric one. This reflects the loss of one of the motion integrals, the 
angular momentum. 
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Figure 4: Three-dimensional orbit for (a) axisymmetric potential and (b) 
triaxial potential, (c) Radial coordinate R as function of time t. (d) Co- 
ordinate z as function of time i. Solid curves represent motion for triaxial 
potential and dashed curves the axisymmetric case. The parameters are the 
same as in figure El Initial conditions: R = 3.0, R = 0, (p = 0, z = 0.5R(p 
and energy E = —0.15. 
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4 Discussion 



By applying a Miyamoto and Nagai transformation on the three Carte- 
sian coordinates of the monopole potential, we obtained a triaxial potential- 
density pair whose mass density distribution is everywhere non-negative 
and free from singularities; also, it is box-shaped with respect to the z 
axis. A triaxial version of one of Satoh's axisymmetric models also yields a 
mass-density distribution with similar characteristics. We believe that these 
simple analytical models may be useful for disk galaxies having box-shaped 
bulges. Some numerically calculated orbits for the Miyamoto and Nagai- like 
triaxial potential were also presented. 

D. Vogt thanks FAPESP for financial support. P. S. Letelier thanks 
CNPq and FAPESP for financial support. This research has made use of 
NASA's Astrophysics Data System. 
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